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Abstract 

We discuss a modified form of gravity implying that the action 
contains a power a of the scalar curvature. Coupling with the cosmic 
fluid is assumed. As equation of state for the fluid, we take the sim- 
plest version where the pressure is proportional to the density. Based 
upon a natural ansatz for the time variation of the scale factor, we 
show that the equations of motion are satisfied for general a. Also 
the condition of conservation of energy and momentum is satisfied. 
Moreover, we investigate the case where the fluid is allowed to possess 
a bulk viscosity, and find the noteworthy fact that consistency of the 
formalism requires the bulk viscosity to be proportional to the power 
(2a — 1) of the scalar expansion. In Einstein's gravity, where a = 1, 
this means that the bulk viscosity is proportional to the scalar ex- 
pansion. This mathematical result is of physical interest; as discussed 
recently by the present authors, there exists in principle a viscosity- 
driven transition of the fluid from the quintessence region into the 
phantom region, implying a future Big Rip singularity. 



1 Introduction 

Recent years have witnessed a considerable interest in modified versions of 
Einstein's gravity, the motivation for this is coming largely from the evi- 
dence that has accumulated about the ongoing accelerated expansion of the 
universe. This evidence is seen in the high redshifts of type la supernovae PP , 
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as well as in the anisotropy power spectrum of the cosmic microwave back- 
ground |2]. These observations lead naturally to the idea of a dark energy, 
believed to be 73% of the total energy content of the universe. Thermody- 
namically, the dark energy concept implies rather bizarre consequences; thus 
the equation of state parameter w of the cosmic fluid may be less than -1, 
implying a singularity in the future called a Big Rip , and it leads also to 
negative entropies In the papers listed in Ref. J3], actually two different 
types of Big Rip were investigated. In the recent paper of Nojiri et al. [S], 
a classification of the known types of Big Rip were made and two new types 
of Big Rip were discovered, so that we know now four types in all. 

The intention of introducing modified gravity models is that one may 
obtain a gravitational alternative to the conventional description of dark 
energy. The possibility may look attractive as the presence of dark energy 
may thereby be a consequence of the expansion of the universe IH]. Various 
aspects of the modified gravity idea have recently been discussed in Refs. jTl 
El m CD] (for an application to the de Sitter space, see Ref. [H].) In the 
present paper we shall study the following modified gravity model: 



Here = SnG, fo and a are constants {a may in principle be negative), and 
Lm is the matter Lagrangian. This kind of system was recently studied by 
Abdalla et al. [H]. Note that the dimension of fo is cm^*^""^). 

It is interesting to note that the string/M-theory effective action may 
produce terms with negative powers of R Now, modified gravity with 
an is in itself not stable ^2] , but the inclusion of a term in the action 
[7j, or the inclusion of quantum effects ^3], makes the theory stable. As we 
are interested in the late time dynamics of our theory we may drop the R^ 
terms and work only with terms that are dominant at late times. 

It turns out that even this apparently simple generalization of Einstein's 
theory is not quite trivial. We will focus attention on the following two 
points: 

• How the equations of motion as derived from the action (^J) are compat- 
ible with a scale factor a varying with time according to the ansatz of Eq. (jUj) 
below. The fluid is in the first part of our work taken to be non-viscous, and 
we adopt for simplicity the conventional equation of state 




(1) 



p= (j-l)p = wp. 



(2) 
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It turns out that the equations of motion are actually satisfied, for arbitrary 
values of a. Moreover, the conservation equation for energy and momen- 
tum are satisfied also. This is reassuring, as regards the adaptability of the 
conventional formalism to the modified gravity theory. 

• Our second point is to introduce a bulk viscosity C in the fiuid. It 
turns out that in order to obtain a consistent description, one arrives in a 
natural way at a time dependent bulk viscosity model; see Eq. fOUj) below. 
In particular, in the case of Einstein's gravity [a = 1), the bulk viscosity 
becomes simply proportional to the scalar expansion. An interesting fact 
is that this is precisely the relationship recently found in Ref. to be 
compatible with a viscosity-driven transfer of the fiuid from the quintessence 
region {w > —1) into the phantom region {w < —1). 

2 Solutions of the Equations of Motion. Non- 
Viscous Case 

We start from the spatially fiat FRW metric. 



and set the cosmological constant A equal to zero. The Hubble parameter 
is H = a/a, where an overdot means differentiation with respect to t. The 
scalar expansion is 6* = f/^;^ = 3H, where is the four-velocity of the fiuid. 

We shall be primarily interested in the development of the late universe, 
from t = to onwards. For simplicity we choose the origin of time such that 
to = 0. The corresponding energy density is po; this quantity is assumed to 
be known. Similarly we assume the initial scalar expansion 6q to be known. 
In Einstein's gravity [111 HZ] 



In the recent modified gravity model of Abdalla et al. the exponent 
2/37 in Eq. is replaced by 20/37. It becomes thus natural to assume the 
following time variation of the scale factor: 



(3) 




(4) 




(5) 



a = ao(l+/3t)'"/='^ 



(6) 
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in the present case. We shall adopt this form in the following. The parameter 
(3 has to be determined from the equations of motion. Taking the (00)- 
component of these equations, it is in principle possible to relate (3 to the 
initial energy density po aX t = 0. Instead of Eq. we may equally well 
write 

« = ao + , (7) 

^0 = 3ifo being the scalar expansion at t = 0. The explicit time dependence 
of 9 thus follows at once, 

^^ = ^^o(l + ^^ot) (8) 

We see that (3 and 9q are equivalent parameters to be determined from the 
equations of motion as they are related through a constant factor, 

/J=^^o. (9) 

In Einstein's gravity, 

/9 = ^7^o, « = 1. (10) 

In this case we find from Riemann's energy equation that 9 is related to p in 
a very simple way, 

9"^ = 3k^p; (11) 

cf., for instance, Refs. IIHIIHI- Equation (fTTj) holds for all values of t. 

From the action ([Q) we obtain in the general case the equations of motion(n] 

+afogf,uV^R''-' = (12) 

where T^i, is the energy-momentum tensor corresponding to (there is a 
printing error before Eq. (3.2) in Ref. [U]). The notation is chosen such that 
the values a = 1, /o = 1 correspond to Einstein's gravity. 
Let us put /i = z/ = 0. In coordinate basis we have 

Roo = , R = q[- + -), 13 



VlR""-^ = (a - l)9o(i?"~'i?), (14) 
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2 DQ-l 



a — 1 



do{a^R"-'R). 



(15) 



Then, since Too = P, the (OO)-component of the equations of motion becomes 



J/oi?" - afo-R''-' - foa{a - l)do{R''-'R) 

Zi Qj 



-afo 



a 



Now taking into account our ansatz (jH)) for the scale factor we calculate 



R 



— 1 



3a 1 ^ (2/j)^"(2a/37-l) a 
a 2(l + /3t)2° 7 

(2/3)2"(a-l)(2a- 1) 



2(1 + /3t)2° 



7 V37 
7 V37 



n o-l 



a-1 



2 oa-l 



-(2/5)'°(a - l)(2a - 2a/7 - 1) 



2(1 + /3t)2° 
whereby we arrive at the equation 



7 V37 , 



a-l 



(16) 

(17) 

(18) 
(19) 
(20) 



/pgg" (7/a)"(4a/37 - 1) 
2(1 + /5t)2° 



a-l 



{(2-«)|^-(a-l)(2a-l)} = «:V, (21) 
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reinstalling ^o by using Eq. ©• As is seen from Eq. (j2H) . the energy density 
varies with time as 



1+I3t) 

If the parameter set {/q, a, po} is given, the value of the initial scalar expan- 
sion 9q can in principle be determined from Eq. (|2ip. setting t = 0. Equation 
(|2ip agrees with Ref. P| (the notation, and the initial assumption (jUI) for the 
scale factor, are different). 

The case of Einstein's gravity, a = 1, /o = 1, yields 



el = 3kVo, 



(23) 



which is in accordance with Eq. (fTT|) when t = 0. It is seen that 9"^ and p are 
both varying with time as (1 + (3t)~'^. 
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Next, we may for general a put fi = u = r in Eq. (|12p. to calculate 
the pressure component of the equations of motion. We have, in coordinate 
basis, Tfj- 

= P9rr = pa ■ We abstain, however, from giving the details here. 
The most important point in our context is that the pressure is from this 
equation found to vary with time in the same manner as the density p, 

Thus, if the equation of state Q holds at t = 0, it holds for all later values 
of t also, with the same constant values for 7 and w. This is important for 
the consistency and applicability of the modified gravity theory. 



2.1 A remark on energy-momentum conservation 

So far, we have not taken into account the condition coming from conser- 
vation of energy and momentum. A characteristic property shows up if one 
takes the covariant divergence on both sides of Eq. ()12j) : terms on the left 
hand side of the equation combine to yield zero. This is shown explicitly in 
a recent calculation of Koivisto Thus we obtain the same conservation 
equation for energy-momentum as in ordinary relativity: 



which in turn implies 



V^T^, = 0, (25) 

(26) 



dp P + P d ^ 



dt a? dt 

Inserting the forms ()22p and for p and p we see that Eq. ()26|) is satisfied 
identically. Again, this is a point supporting the consistency of the modified 
gravity. 

We may moreover note that the relation 



Po VOo 



(27) 



known from ordinary relativity, is satisfied also here, for arbitrary values of 
a. 
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3 Introduction of Viscosity 



Assume now that the fluid is viscous. Because of the assumed spherical 
symmetry the shear viscosity plays no role, and only the bulk viscosity ( has 
to be considered. On thermodynamical grounds, in conventional physics ( 
has to be positive; this being a consequence of the positive entropy change 
in irreversible processes. We shall assume C > ^ilso here. In general, the 
presence of ( does not have any influence upon the (OO)-component of the 
equations of motion. The only change in the formalism because of viscosity 
is that the thermodynamical pressure p becomes replaced with the effective 
pressure p, defined as 

p = p-C9 (28) 

(cf., for instance, Refs. fHl E!)- This follows from the expression for the 
stress tensor for a viscous fluid. It can be easily shown that the (rr)- 
component of Eq. (fT^ (the pressure component of the equations of motion) 
becomes satisfied for all t if C is taken to be time dependent, ( oc 6*^°"^, or 

^ " (l + /3t)2"-i' ^^^^ 
Namely, the product (6 then varies with t in the same manner as p, 

cf. Eq. (1211). We do not need here to involve the pressure equation explicitly, 
but can confine ourselves to the energy equation (j2H), making use of the 
equation of state In the presence of viscosity the right hand side of 
Eq. (j?T|l thereby becomes changed to 

K^p = -^p ^ -^{p - ce). (31) 

7 — 1 7 — 1 

Thus the generalized equation from which 6o can be determined in the viscous 
case, becomes 

I |(2-«)^-(a-l)(2a-l)|c 

2 / _ Co ^0 
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Let us here consider two special cases. As usual, the most important case 
is that of Einstein's gravity, a = 1, /o = 1, whereby Eq. (j22I) reduces to a 
quadratic equation in ^o- We get 



-3^X0 + a + 12(7 - lypo/^' 
^° = 2i^) • ^ ^ 

The case a = 1 turns out to have a bearing on the transition of the fluid 
through the barrier w = 1 into the phantom region. In the notation of 
Ref. we can write in this case 

C = re, (34) 

r being a positive constant. As found in Ref. [12], if the magnitude of the 
bulk viscosity is so large that the condition 

3kV > 7 (35) 

is satisfied, then the fluid will be driven into the phantom region {w < —1) 
and the Big Rip singularity, even it started out in the quintessence region 
{w > —1). This result follows directly from Friedmann's equations for a 
viscous fluid, and shows up mathematically as a singularity in the density p 
occurring at a finite time t. 

As our second example we will consider the case a = 1/2, meaning that 
the gravitational action involves the square root of the scalar curvature. 
Mathematically, Eq. (jH^ then becomes quite simple and yields the solution 

e, = ^ — --. (36) 

Joy 16(1-37/2) 7-1 

The initial scalar expansion becomes thus proportional to the initial density. 
The same proportionality holds true for the time dependent quantities 9 and 
p; cf. Eqs. © and ((221). 



4 Summary 

We started out with a spatially flat FRW space, adopting a modified gravity 
action in the form of Eq. (^Q); this form being recently studied in Ref. 0. The 
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corresponding equations of motion are given in Eq. (jl2p . We made the basic 
ansatz (jUI) for the time variation of the scale factor a. This ansatz turned 
out to function well throughout, thus showing in general the flexibility and 
usefulness of the modified gravity formalism. We considered also the energy- 
momentum conservation, with similar satisfactory results. 

As far as viscosity is concerned, the formalism admits the presence of a 
bulk viscosity C, without any problems if C, is taken to be proportional to 
the scalar expansion raised to the power (2a — 1); cf. Eq. (j^ . Then 
the product C, 9 varies with time in the same manner as the pressure p, and 
the equations of motion are satisfied. This behavior could hardly have been 
seen beforehand, without a detailed consideration of the formalism. In the 
special case of Einstein's gravity (a = 1) this result fits nicely into the viscous 
cosmology description recently given in Ref. there exists in principle a 
viscosity-driven transition of the fluid from the quintessence region into the 
phantom region, implying a future Big Rip singularity. 

Finally we mention that our inclusion of bulk viscosity may be used in 
more general modified gravities, for instance, those introduced by Allemandi 
et al. dHl. 
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